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ABSTRACT 

Strong gravitational lensing is a powerful tool for probing the matter distribution 
in the cores of massive dark matter haloes. Recent and ongoing analyses of galaxy 
cluster surveys (MACS, CFHTLS, SDSS, SGAS, CLASH, LoCuSS) have adressed the 
question of the nature of the dark matter distribution in clusters. N-body simulations 
of cold dark-matter haloes consistently find that haloes should be characterized by a 
concentration-mass relation that decreases monotonically with halo mass, and popu- 
lated by a large amount of substructures, representing the cores of accreted progenitor 
halos. It is important for our understanding of dark matter to test these predictions. 
We present MOKA, a new algorithm for simulating the gravitational lensing signal from 
cluster-sized haloes. It implements the most recent results from numerical simulations 
to create realistic cluster-scale lenses with properties independent of numerical resolu- 
tion. We perform systematic studies of the strong lensing cross section as a function of 
halo structures. We find that the strong lensing cross sections depend most strongly 
on the concentration and on the inner slope of the density profile of a halo, followed 
in order of importance by halo triaxiality and the presence of a bright central galaxy. 

Key words: galaxies: halos - cosmology: theory - dark matter - methods: analytical 
- gravitational lensing: strong 



1 INTRODUCTION 

Galaxy clusters are an important probe for dark-matter 
properties. According to the standard scenario of structure 
formation, galaxy clusters as a population are still in their 
formation process. Since gas cooling cannot substantially 
compress dark matter haloes, their density profiles are dom- 
inated by dark matter. 

Recent analyses of strong and weak lensing by galaxy 
clusters have found good consistence with an NFW-like pro- 
file whose projected mass density distribution continuously 
flattens towards the centre, as expected from CDM dom- 
inated haloes. However, for an increasing number of clus- 
ters Oguri et al. (2005, 2009); Umetsu et al. (2011) have 
found that the mass-concentration relation seems to lie sub- 
stantially above the relation predicted by CDM simulations. 
Strong lensing by clusters is sensitive to their internal struc- 



ture - mass distribution within the Einstein radius , which 
includes (i) presence of substructure, (ii) asymmetry in their 
gravitational potential well, (iii) ellipticity, (iv) presence of 
a massive and bright central galaxy (v) inner slope of the 
dark matter density profile. 

Apart from weak and strong lensing, the halo density 
profile can also be constrained from the velocity dispersion 
profile of the central galaxy (Sand et al. 2004; Newman et al. 
2009) and X-ray emission from the hot intra-cluster gas. 
While gravitational lensing does not rely on any equilibrium 
or dynamical assumptions, methods based on galaxy or gas 
dynamics do, potentially biasing mass and concentration es- 
timates. 

In this paper, we shall quantify the importance of the 
structural parameters of dark matter haloes for strong lens- 
ing signal, using our new and fast algorithm MOKA (Matter 
density distribution Kode for gravitationAl lenses) to create 
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1 For an axially symmetric lens is defined as the radius of a circle 
enclosing a mean convergence of 1. 
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realistic maps of substructured triaxial dark matter haloes, 
which are in perfect agreement with the measurements done 
on galaxy clusters extracted from numerical simulations. 

In Section 2 we present the halo properties on which our 
algorithm relies. The halo lensing properties are presented 
in Section 2.2 and the strong lensing signal dependence on 
the halo modelling is discussed in Sections 3. Discussion and 
conclusion are presented in Section 4. 

We adopt a ACDM model with tt m = 0.3, Q A = 0.7, 
h = 0.7 and as = 0.9, consistently with the Mare Nostrum 
Universe that we compare some of our results to. 



2 CONSTRUCTION OF REALISTIC LENSES 

Strong gravitational lensing depends on the matter distri- 
bution near halo centres where limited numerical resolu- 
tion may affect the particle distribution. Here, we present a 
new algorithm, MOKA, which analytically creates surface mass 
density distributions of triaxial and substructured haloes in- 
dependent of numerical resolution. MOKA is publicly available 
at the following url http:/ /cgiocoli. wordpress.com/research- 
interests/moka. The idea behind this new algorithm is to 
construct realistic lenses starting from a set of ingredients 
taken from state of the art numerical simulations. It not 
only accounts for the smooth dark halo and stellar compo- 
nents, as earlier studies do (van de Ven et al. 2009; Man- 
delbaum et al. 2009), but also the presence of substructures 
perturbing the regular matter distribution. The procedure 
is discussed in the following subsections. 



2.1 Host Halo Contents 

Halo properties are described by the extended halo model 
Giocoli et al. (2010), developed for the reconstruction of the 
non-linear dark matter power spectrum. We briefly list its 
ingredients here. 



2.1.1 Properties of the dark matter component 

The virial mass of a halo is defined as 
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where p c represents the critical density of the Universe, fio = 
fim(O) the matter density parameter at the present time 
and A„i r is the virial overdensity (Eke et al. 1996; Bryan 
& Norman 1998), R V i r symbolizes the virial radius of the 
halo which defines the dinstance from the halo centre that 
encloses the desired density contrast. In what follows, we 
summarize the most recent numerical results built into MOKA. 

• As said above, the dark matter density distribution in 
isolated haloes is well described by the NFW (Navarro et al. 
1996) profile 
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where r 3 is the scale radius, defining the concentration 
c vir = R vir /r s , and p s is the dark matter density at the 
scale radius, 
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Combining the preceding three equations, we can explicitly 
write the NFW profile as a function of c V i r and M V i r , 
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We use this density profile to model the dark matter distri- 
bution in the clusters produced by MOKA. 

• The halo concentration is a decreasing function of the 
host halo mass. This numerical result is explained in terms 
of hierarchical clustering in CDM and the different halo- 
formation histories. Already Bond et al. (1991) and Lacey 
& Cole (1993), following the extended-Press & Schechter 
(1974) theory, have found that the collapse time of dark mat- 
ter haloes depends on the halo mass and that their assembly 
history is hierarchical: small systems collapse at higher red- 
shifts than larger ones (Sheth & Tormen 2004a; Giocoli et al. 
2007). This trend is reflected in the mass-concentration rela- 
tion: at a given redshift smaller haloes are more concentrated 
than larger ones. Different fitting functions for numerical 
mass-concentration relations have been given (Bullock et al. 
2001; Neto et al. 2007; Duffy et al. 2008; Gao et al. 2008). In 
this work, we use the relation proposed by Zhao et al. (2009) 
which links the concentration of a given halo with the time 
(to. 04) at which its main progenitor assembles 4 percent of 
its mass, 
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The model by Zhao et al. (2009) fits well numerical sim- 
ulations even with different cosmologies. It seems to be of 
reasonably general validity. 

Due to different assembly histories, haloes with same 
mass at the same redshift may have different concentrations 
(Navarro et al. 1996; Jing 2000; Wechsler et al. 2002; Zhao 
et al. 2003a, b). At fixed host halo mass, the distribution 
in concentration is well fitted by a log normal distribution 
function with a variance oi nc between 0.1 and 0.25 (Jing 
2000; Dolag et al. 2004; Sheth & Tormen 2004b; Neto et al. 
2007). 

• Halos are generally not spherical but triaxial due to 
their tidal interaction with the surrounding density field dur- 
ing their collapse (Sheth & Tormen 1999, 2002). The corre- 
lation of the halo shape with the surrounding environment 
is expressed as a function of the matter density parameter 
and of the typical collapse mass. Jing & Suto (2002) have 
performed a statistical study of halo shapes extracted from 
numerical simulations, deriving that, if a, b and c represent 
the minor, median and major axes, respectively, empirical 
relations for a/c and a/b are given by the distribution 
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for A = (a/c)(M„ il ./M t (z ! )) a07n(Z!) with a x 
the conditional probability for the axis ratios 
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0.5 if a/c < 0.5 else r m i n = a/c. As usual, 
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M*(zi) is the non-linear mass at redshift z\. Simulations 
with gas cooling have produced haloes that tend to be more 
spherical than those in pure DM simulation (Kazantzidis 
et al. 2004). The host halo shape then depends also on the 
morphology of the stellar mass component. In this paper 
we use the prescriptions by Jing & Suto (2002), which are 
based on dark matter only simulations. However, MOKA is a 
very flexible tool, and the input distributions of axis ratios 
can be easily modified. 

• Haloes are not smooth, but characterized by a large 
number of substructures, which may or not host satellite 
galaxies (Moore et al. 1999; Springel et al. 2001b; De Lu- 
cia et al. 2004). These substructures are cores of progenitor 
haloes accreted along the merger tree that were not com- 
pletely disrupted by tidal stripping (Springel et al. 2001b; 
Gao et al. 2004; De Lucia et al. 2004; van den Bosch et al. 
2005; Giocoli et al. 2008; Shaw et al. 2007). Because of dif- 
ferent assembly histories and time scales for subhalo mass 
loss, more massive haloes retain more substructures than 
less massive haloes at a given redshift. Likewise, haloes with 
a lower concentration (and thus with a lower formation red- 
shift) are on average more substructured than haloes with 
larger concentration. Giocoli et al. (2010a) fitted the subhalo 
mass function by 
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where A = 9.33 x 10" 4 , /3 = 12.2715, a = -0.9 and c is the 
mean concentration of a halo with mass M v i r at redshift 
zi. We adopt this model because it incorporates the depen- 
dences of the subhalo mass function on the mass M vir , the 
concentration c V i r and the redshift zi of the host halo. To 
populate our haloes with substructures with mass m, we 
will randomly sample the distribution down to a minimal 
subhalo mass m m i n . 

Subhalo density profiles are modified by tidal stripping 
due to close interactions with the main halo smooth com- 
ponent and to close encounters with other clumps, gravi- 
tational heating, and dynamical friction. Such events can 
cause the subhaloes to lose mass, and may eventually result 
in their complete disruption (Hayashi et al. 2003; Choi et al. 
2007). The remaining self-bound subhaloes will have density 
profiles different from the NFW shape in that they are trun- 
cated at the tidal radius (Bullock et al. 2001). In order to 
take the truncation into account, we model the dark matter 
density distribution in subhaloes using truncated singular 
isothermal spheres (Keeton 2003), 

f < R 

Psub{r) = { 2 ^Gr 2 r Hsub > 
[0 r > R sub 

with velocity dispersion u v , and R S ub defined as: 
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Based on the temperature definition in the spherical collapse 
model, we can write 
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where k is Boltzmann's constant, and M15 is the mass in 
units of 10 15 Mq//i. For temperatures much above 1 keV, 
Eq. (9) follows the observed relation between X-ray temper- 
ature and velocity dispersion by Wu et al. (1998), 
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while for T <gil keV it reduces to the relation 
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that can be obtained combining Eqs. (8) and (1) identifying 

Rsub — Rvir- 

Let us define R t as the radius at which the subhalo mean 
density is of the order of the mean density of the main halo 
within r. Following Tormen et al. (1998) we can write it as: 
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The velocity dispersion a v is related to the subhalo temper- 



where m S ub an d r represent the subhalo mass and its dis- 
tance from the host halo centre, and M„; r (r) the host halo 
mass profile. Truncating subhaloes at R S ub does not cre- 
ate discontinuities in the convergence map because Rt(r) < 
Rsub- It preserves the total subhalo mass fraction in hosts as 
found in numerical simulations (Gao et al. 2004; De Lucia 
et al. 2004; Giocoli et al. 2010a). If the host halo matter den- 
sity distribution is described by an NFW- profile, the subhalo 
tidal radius as a function of the distance from the host halo 
centre can be analytically estimated. In Figure 1 we show 
the ratio between the subhalo tidal radius and the radius 
which encloses the subhalo mass as a function of the dis- 
tance from the host halo centre. From the figure we notice 
that for r < 0.75R V i r Rt <C Rsub while when the subhalo is 
located near at the virial radius the ratio tends to unity. 

The SIS profile profile well represents galaxy density pro- 
files on scales relevant for strong lensing. Previously, differ- 
ent authors have used this model to characterize the lensing 
signal by substructures after stripping (Metcalf & Madau 
2001). Nonetheless, additional subhalo density profiles will 
be implemented in MOKA to allow truncating the subhalo pro- 
file more smoothly. 

• The spatial distribution of subhaloes tends to follow the 
smooth dark matter distribution of the host halo. However, 
clumps near the host halo center are easily destroyed, due 
to their tidal interaction with the main halo. Thus, the spa- 
tial subhalo distribution is less concentrated than the NFW 
profile of the host halo. Gao et al. (2004) have studied this 
using a cosmological numerical simulation and found that 
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Figure 1. Ratio between the subhalo tidal radius and the one 
which encloses its mass as a function of the distance from the 
host halo centre. In the figure we show the case of three subhaloes 
located in a host halo with mass 10 15 M e /h at redshift z = 0.25 
with concentration c v i r = 4. 



the cumulative spatial density distribution of clumps in host 
haloes is well described by 



n(< as) (1 + a'c vir )2 



(15) 



Ntot (1 + a'c vlr x 2 ) 

where x is the distance to the host halo centre in unit of 
the virial radius, Ntot is the total number of subhaloes in 
the host, a' = 0.244 and /3' = 2.75. In the Figure 2 we 
show previouse equation for different value of the host halo 
concentration, as indicated in the label. To place sublumps 
in the host, we sample this distribution and randomly assign 
position angles 6 and on a sphere. A similar equation has 
been used to model the satellite galaxy density distribution 
by van den Bosch et al. (2004). 



2.1.2 Dissipative 
dark matter 



nic component and its effects on 



Strong lensing is sensitive to the matter distribution inside 
the centre of galaxy clusters (r ~ 100 kpc). On such scales, 
the density of the baryonic component becomes comparable 
to that of the dark matter. Meneghetti et al. (2003) have 
shown that the brightest central galaxy (BCG) on strong 
lensing signal is generally moderately important for strong 
lensing, but potentially decisive for those haloes which are 
otherwise marginally supercritical strong lenses. 

To populate haloes with a central galaxy of a cer- 
tain stellar mass, we use the Halo Occupation Distribution 
(HOD) technique. HOD assumes that the stellar mass of a 
galaxy is tightly correlated with the depth of the potential 
well of the halo within which it formed, thus 
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Figure 2. Radial subhalo density distribution model. The differ- 
ent curves show equation (15) for different values of the host halo 
concentration c v i r , as in the label. 
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Figure 3. Density profile of a halo at 2; = 0.25 with M v i r 

■ The dotted curve shows the NFW profile, the dashed 
curve the Henquist profile for the stellar component of the BCG 
at the halo centre with stellar mass M a tar = 5 X 10 12 Mq /h. The 
dissipative baryon component modifies the dark-matter profile by 
adiabatic contraction as shown by the long-dashed line. The solid 
line shows the final profile. 



as estimated by Wang et al. (2006), who modelled this re- 
lation after the semi-analytic galaxy catalogue of the Mil- 
lennium Simulation (Croton et al. 2006). In this relation, 
we include a Gaussian scatter in M s tar at a given host halo 
mass with (JM ataT = 0.148. For the central galaxy, we set the 



parameters a" = 0.39, = 1.96, and log (M, 



tarfi) 



10.35. 



For the stellar component of the BCG residing in the 
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halo centre, we adopt the Hernquist (1990) profile 



Pstar(r) 
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It has a scale radius r g related to the half-mass (or effective) 
radius R e by r g = 0.551-R e , as done by Keeton (2001) we 
define the effective radius to be R e — 0.03i?„i r . The scale 
density p g can be estimated by the definition of the total 
mass of a Hernquist model, 



M s tar 
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The presence of a dissipative baryonic component influ- 
ences the dark matter distribution near the host halo centre. 
Blumenthal et al. (1986) described the adiabatic contraction 
analytically, finding good agreement with numerical simula- 
tions. The initial and final density profiles - characterized 
by an initial radius and a final radius r/, when a central 
galaxy is present - are related by 



' [M star (r) + Mdmj] = nM DM ,i{ r i) , (19) 



where 
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Figure 4. Convergence profile of the same cluster, assuming 
sources at rcdshift z s =2. 



M DM j = M DM ,i (1 - fcooi) , (20) 

and f C ooi is the baryon fraction in the halo that cools to 
form the central galaxy. To solve the adiabatic-contraction 
equation, we need to derive r from equation (19). With a 
Hernquist model, 

fcooir 3 + (r + r g ) 2 [(1 - f coo i)r - n] m l {r j ) = ; (21) 

This equation has a single relevant real root. We recall that 
rrii in the previous equation defines the initial mass pro- 
file normalized by the virial mass. Figure 3 shows the den- 
sity profile of a halo with virial mass M v i r = 10 M@/h 
and c V ir = 4 populated with a galaxy with stellar mass 
M Bt ar = 5 x 1O 12 M0 /h. The dotted and short-dashed curves 
refer to the initial dark matter and stellar density in the 
halo, respectively. The presence of a dissipative baryonic 
component contracts the dark matter distribution. Solving 
Eq. (21), the dark matter density profile changes to the 
long-dashed curve. The solid curve shows the total density 
distribution obtained by summing the short and the long- 
dashed curves. We recall that for DM-only realizations, M v i r 
is given by the sum in smooth plus clumpy components, 



while M V ir = 
present. 



M a 



+ Ei=T mi + M a tar if a BCG is 



2.2 Halo Lensing Properties 

In this section we will describe how we calculate lensing 
properties starting from the 3D matter density of all com- 
ponents characterising the haloes. For each component, we 
project the density on a plane perpendicular to the line of 
sight, 



T,(x,y) 
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p(x,y,z)dz. 
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where r = \J x 2 + y 2 + z 2 with the coordinate origin put 
into the host halo centre, and defining £ = ^Jx 2 + y 2 . 

The projected mass density of the NFW, Hernquist and 



SIS profiles can be given analytically (Bartelmann 1996 
Bartelmann & Schneider 2001; Keeton 2001): 
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where £ = £/r s , rj = £/r g , and the two functions F(£) and 
G{r)) are 
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C> l, 

C<1, 
C = i; 



V > 1, 



tan 1 \Jr\ 2 — 
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We recall that the subhalo density profile is truncated at the 
radius enclosing its total (bound) mass. The convergence is 
the appropriately scaled surface mass density y), 



n(x,y) 



where 



D a 



4nG DiDu 



(26) 



(27) 



is the critical surface mass density, containing the angular- 
diametre distance Di, D s and Di s from the observer to the 
lens, to the source, and from the lens to the source, respec- 
tively. 

Figure 4 shows the convergence profiles for the halo 
components presented in Figure 3, with sources assumed 

3jtj Zg — — 2 . 
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Figure 5. Convergence maps of cluster-sized haloes created with our MOKA code at redshift z\ = 0.25, assuming sources at the single 
redshift z a = 2. The solid curves indicate the tangential and radial critical curves. All haloes have a virial mass of 10 15 Mq//i and a 
minimum subhalo mass of 10 10 Mq/Zi. 



For each dark-matter halo, the total convergence map 
is the sum of all contribution maps, 



k(x, y) = kdm{x, y) + 



(28) 



Kstar(x,y) + ^ K su b,i (x — x c , y — y c ) , 



where x c and y c represent the center of mass of the i-th sub- 
structure. To introduce ellipticity into our model, we draw 
random numbers from the Jing & Suto (2002) distributions 
for the axial ratios a/b and a/c, requiring abc — 1. Once 
the axial ratios are known, we deform the convergence map 
accordingly. To randomly orient the halo, we choose a ran- 
dom point on a sphere identified by its azimuthal and eleva- 
tion angles and rotate the halo ellipsoid by these angles. We 
assign the same projected ellipticity to the smooth compo- 
nent, to the stellar density and to the subhalo spatial den- 
sity distributions. For the dark matter density distribution 
in subhaloes we assume spherically symmetric models, this 
because the subhalo typical scale is much smaller than the 
virial radius of the host system in which they are located. 
Having elliptical subhaloes will impact slightly on the total 
strong lensing cross-section. We have chosen to assign the 
same ellipticity to the BCG and to the dark matter den- 
sity distribution in the host following the results of Fasano 
et al. (2010) who found that the shape of the BCG tends 
to reflect that of the associated dark matter halo. Figure 
5 show the convergence maps of six galaxy clusters gener- 
ated with our algorithm. All haloes are located at redshift 
z — 0.25 and possess a virial mass equal to 10 15 Mq / h. The 
source redshift is z s = 2. The subhalo mass resolution is 
10 10 M Q /h. This ensures a substructure mass fraction com- 
patible with Richard et al. (2010), and consistent with the 



fact that systems with a lower substructure fraction tend 
to form at higher redshift, and are thus more concentrated 
(Smith & Taylor 2008). 

From the convergence map, the effective potential and 
the scaled deflection angle can be obtained by means of the 
equations 



9(x,y) = - f K(f)m|f-|V 2 £' 



and 



a(x,y) 



~ [ K 



^4^4dY 
if -ei 



(29) 



(30) 



A source will be seen at the angular position 9 related to its 
intrinsic angular position by the lens equation 



(31) 



Derivatives of the lensing potential are denoted by sub- 
scripts, 

d 2 $(x,y) 



(32) 



where £ j = x when i = 1 and & = y when i = 2. We 
introduce the pseudo- vector field of the shear 7 = (71,72) 
by its components 



!i(x,y) = -($11 -$22) 

72(2;,?/) = $12 = $21- 



(33) 
(34) 



Light bundles are deflected differentially and thus dis- 
torted as described by the Jacobian matrix 



1 'i'.,)- 



(35) 
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with the eigenvalues 

A t = 1 — k — 7, (36) 
X r = 1-K + f. (37) 

The conditions At = and A r = define the location of the 
tangential and radial critical curves in the lens plane, where 
the magnification is formally infinite. Mapping the critical 
curve back into the source plane by the lens equation gives 
the tangential and radial caustics. 

Table 1 summarizes the list of parameters controlling 
MOKA. 



3 LENS STRUCTURAL PROPERTIES AND 
STRONG LENSING SIGNALS 

In this section we discuss the strong lensing efficiency of 
clusters produced by MOKA. Most of the discussion is focussed 
on the lensing cross section a. This is defined as the region 
on the source plane from where the sources are mapped into 
images with a certain length-to- width ration l/w. 

For each cluster we create a high-resolution deflection 
angle map of 2048 x 2048 pixels centred on the cluster cen- 
tre, with side length equal to the virial radius. From this 
we numerically estimate the lensing cross section by means 
of ray-tracing methods (Meneghetti et al. 2005a, b; Fedeli 
et al. 2006; Meneghetti et al. 2011). A typical calculation 
takes about one minute on a 3.06 GHz single-core proces- 
sor. Starting from the observer, bundles of light rays are 
traced back to the source plane. Using an adaptive grid, 
this is populated by elliptical sources of a fixed equivalent 
radius of 0.5". The number of highly magnified images is 
increased by refining the spatial resolution on the source 
plane near caustics. Analysing the images individually and 
measuring their length- width ratio l/w, we can define the 
lensing cross section for giant arcs &i/ w as the region on the 
source plane from where the sources are mapped into im- 
ages with a certain l/w. Giant arcs are commonly defined as 
distorted images with a length-to- width ratio l/w ^ 7.5 or 
l/w 10. We shall discuss our results for the three different 
values l/w = 5, 7.5 and 10. 

The integral of the strong lensing cross section of a 
galaxy cluster for different z s , weighted by a source red- 
shift distribution, allows to quantify the number of gravita- 
tional arcs expected to be seen. The strong lensing analysis 
of galaxy cluster size haloes from a ACDM numerical simu- 
lations by Bartelmann et al. (1998) has revealed that those 
haloes produce an order of magnitude fewer gravitational 
arcs than observed. Studying how strong lensing cross sec- 
tions change with structural halo properties may help to 
understand possible current limitations of numerical simu- 
lations, without the need of advocating cosmological models 
with dark energy (Bartelmann et al. 2003). 

Recent observational studies of strong lensing clusters 
revealed that observed clusters seem to be more concen- 
trated than simulated clusters. This causes more prominent 
strong lensing features. In the following, we shall quantify 
how much the presence of a BCG and halo triaxiality influ- 
ence the strong lensing cross section for giant arcs. We have 
created a sample of 128 high-resolution convergence maps of 
haloes with mass 10 15 M Q /h and sources located at a fixed 



redshift z s — 2. This study is considered as a reference start- 
ing point to present our code MOKA. At the present time we 
are performing a statistical analysis of haloes in a ACDM 
universe to give an estimate of the strong lensing cross sec- 
tion as a function of redshift and of the halo abundance 
(Boldrin et al. in preparation). 



3.1 BCG 

Using a variety of Gas-dynamical simulations, Puchwein 
et al. (2005) studied the impact of the gas component on 
strong lensing signals, finding that the formation of stars 
mainly in the central region of the cluster tends to in- 
crease the cross section by about 20 — 30%, depending on 
lens redshift. To consider an analogous case, we have gen- 
erated a sample of haloes with the same structural proper- 
ties (Mvir = W 15 M Q /h, zi = 0.25 and z s = 2), with and 
without a central galaxy. We recall that for pure dark mat- 
ter lenses the virial mass takes into account the sum of the 
smooth plus clump components, while in simulations with 
BCG the total M v i r includes also the presence of the stellar 
system. In both cases, the total masses are identical. For 
each halo we have post-processed their deflection angle map 
and estimated the strong lensing cross sections as a function 
of arc length-to- width ratio l/w. The left panel of Figure 6 
shows the median strong lensing cross section as a function 
of l/w for a sample of clusters with (solid) and without BCG 
(dashed curve). The shaded regions enclose 25 — 75% of the 
data. The right panel of the same figure shows for each clus- 
ter the ratio of the strong lensing cross section as function 
of length-to- width ratio with and without a BCG, including 
the estimated median. The Figure confirms that at redshift 
zi =0.25, the value of a for arcs with length-to- width ratio 
l/w > 5 tends to be lower of a 20% when the BCG is not 
included. 



3.2 Triaxiality 

Now we consider the case of spherical and triaxial haloes. 
Triaxiality affects strong lensing through the halo orienta- 
tion. A cluster whose major axis is oriented along the line 
of sight will be a more efficient strong lens than if it is ori- 
ented otherwise. This orientation bias influences also mass 
and concentration estimates (Meneghetti et al. 2010b). To 
quantify how important triaxiality is for the appearance of 
giant arcs, we have created a sample of spherical haloes and 
compared their strong lensing cross sections to those of our 
fiducial sample of triaxial haloes whose mean ellipticity is 
£313 = 0.15 . Figure 7 shows the strong lensing cross section 
as a function of l/w for the same sample of either triaxial or 
spherical galaxy-cluster haloes. Both samples contain BCGs. 
The right panel shows the median ratio of a between the two 
samples: for l/w > 5, the strong- lensing signal tends to de- 
crease by between 20 % and 70 % in spherical compared to 
triaxial haloes. This confirms results obtained with a small 
sample of simulated haloes by Meneghetti et al. (2007b). 



2 We define the 3D ellipticity as: = (c—a)/[2(a+b+c)], where 
a, b and c define the minor, median and major axes respectively. 
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Figure 6. Left panel: median of the strong lensing cross sections a of a sample of simulated haloes with 10 15 Mq//i at 2; = 0.25 with 
sources at z s = 2 as a function of l/w. The solid curves show the median and 25 — 75% percentiles of a for triaxial haloes with a BCG, 
while the dashed curves are the same sample but with only DM components. Right panel: median ratio of the strong lensing cross section 
without and with BCG as function of l/w. The scatter for l/w > 5 is of the order of 0.4. 



Table 1. List of parameters present in MOKA 1.0: http://cgiocoli.wordpress.com/research-interests/moka. 

Halo concentration: drawn from a log-normal distribution around the mean value at fixed mass 1 

Halo virial radius: defined from the spherical collapse model 

Axial ratios: drawn from the Jing & Suto (2002) model 2 

Halo orientation: random picking the point on a sphere 

Subhalo mass function: random sampling the Giocoli et al. (2010a) model 

Subhalo distribution: random sampling the Gao et al. (2004) model, also NFW model has been implemented 

Subhalo velocity dispersion: fixed by the equation (9) 

BCG stellar mass: sampling the Wang et al. (2006) model with a gaussian scatter 

BCG effective radius: fixed to be the 4% of the host halo virial radius (Kceton 2001) 



x M0KA is flexible being able to work with different mass-concentration relation model 
2 cllipticity can be turned off 



Notice also that arcs with small length-to-width ratios (be- 
tween 2 to and 5) are more abundant in spherical haloes 
because of their relatively lower shear. 



each halo Fedeli et al. (2010) estimate the strong lensing 
cross section as mean of the three projections: at the end 
this procedure cuts the wings of their distribution. 



3.3 Scatter 

At a given host halo mass and lensing redshift, the strong 
lensing cross section scatters around its mean. This scat- 
ter is due to both halo triaxiality, presence and distribution 
of substructures and concentrations. Fedeli et al. (2010) fit 
the scatter in y — log(a7.s) for galaxy cluster sized haloes 
with a Gaussian distribution, or an Edgeworth expansion 
around it, including skewness. Figure 8 shows the distribu- 
tion of the strong lensing cross sections for three length-to- 
width ratios l/w = 5, 7.5 and 10. In all cases, Gaussian 
distributions fit the data well with no significant deviation 
from the Edgeworth expansion, which we superpose. The 
distributions have standard deviation of ffiogfo-) slightly in- 
creasing with l/w; we find ci og ( CT ) = 0.56, 0.59 and 0.61 for 
l/w — 5, 7.5 and 10, respectively. These values are approxi- 
mately twice those found by Fedeli et al. (2010). We believe 
that the main reason for the larger standard deviation is how 
the halo boundaries were defined in the simulations and with 
our algorithm: while we include all matter of the main halo 
along the line of sight, Fedeli et al. (2010) have extracted 
boxes of 5 Mpc side length from their simulation. Another 
contribution to this difference is given by the fact that for 



3.4 Concentration 

Haloes of the same mass have different concentrations re- 
flecting their different assembly histories. We show in Ap- 
pendix A how the rms scatter in the strong lensing cross 
sections is related to the concentration scatter and thus to 
the different host halo merger histories. 

Strong lensing depends mainly on the matter density 
in the cluster cores, where the dark-matter density profile 
approaches a logarithmic slope near —1. The scale where 
the logarithmic density slope is —2 defines the scale ra- 
dius r s and the host halo concentration by c V i r — r s /R V i r . 
Higher concentrations cause larger strong lensing cross sec- 
tions. Figure 9 shows this relation for our halo sample at 
zi = 0.25. Filled triangles show the median of the relation, 
with error bars enclosing the quartile of the distribution in 
each (77.5 bin. The solid line shows a least-squares fit to the 
data whose slope is 2.83 ± 0.06. The dashed curve shows 
the predicted relation between the strong lensing cross sec- 
tion and the host halo concentration for spherical smooth 
NFW haloes with mass 10 15 M Q //i. Notice that our simu- 
lated lenses tend to have larger strong lensing cross sections 
at low concentration compared to spherical NFW haloes. 
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Figure 7. Left panel: median of the strong lensing cross sections a of a sample of simulated haloes with 10 15 Mq//i at z\ = 0.25 with 
sources at z s = 2 as function of l/w. The solid line shows the median and the 25 — 75% percentiles of a for triaxial haloes (same as 
Fig. 6), while the dashed line is for spherical haloes. BCGs are present in both samples. Right panel: median of the ratio of the strong 
lensing cross section for the trixial and spherical haloes as function of l/w. Neglecting triaxiality decreases the singal for l/w > 5 by 
20 - 70%. The scatter for l/w > 5 is of the order of 0.4. 
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Figure 8. Distribution of the strong lensing cross section at 
l/w = 5, 7.5 and 10. To fit the distributions, we have considered 
both a Gaussian and an Edgeworth expansion, for y = log(cr), 
which are superposed for corresponding value of cry w . The solid 
line shows a Gaussian with a = 0.3. 

This is mainly due to the BCGs and triaxiality which both 
tend to increase the projected central mass density. 

3.4-1 Statistical Properties 

To summarize, we discuss now the strong lensing efficiency of 
a population of clusters generated with MDKA and we compare 
it to that of a similar sample of numerically simulated cluster 
taken from a cosmological box. 

Meneghetti et al. (2010a) studied the strong lensing by a 
sample of galaxy clusters from the MARE NOSTRUM UNIVERSE. 
In Figure 10, we compare our findings with the median cross 
section of a sample of simulated haloes with masses between 
6 - 7 x W 14 M Q /h at redshift zi = 0.18 (filled circle). Since 
the MARE NOSTRUM simulation does not include radiative pro- 
cesses and star formation, haloes formed therein do not con- 
tain a central massive galaxy. To make them compatible 
with the M0KA haloes, we add a BCG with mean 10 12 M Q /h 
at their center. The solid line shows the median cross sec- 
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Figure 9. Correlation between the strong lensing cross section 
for arcs with length-to-width ratio ^ 7.5 and the host halo con- 
centration. Filled triangles show the median of the correlation and 
the shaded region encloses 25 — 75% of the data. The solid line is 
a least-squares fit to these data points whose slope is 2.83 ± 0.06. 
The dashed curve shows the prediction for spherical and smooth 
NFW haloes (without a central galaxy). In this case the sample of 
haloes is made up of 1024 systems with mass M vir = 10 15 M @ /h 
at redshift z\ = 0.25. 

tion (shaded regions enclose 25 — 90% and 90 — 100% of the 
data) for a sample of M0KA haloes in the same mass range, 
whose abundance has been weighted by the Sheth & Tormen 
(1999) mass function. The three solid curves show the same 
regions for the haloes in the MARE NOSTRUM UNIVERSE. Our 
algorithm well reproduces the median of the strong lens- 
ing cross sections measured in the numerical simulation. It 
doesn't reproduce the scatter of the lensing cross sections 
for the reasons discussed above. 
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Figure 10. Comparison between the median strong lensing cross 
section a as a function of l/w measured in analytically modelled 
haloes (solid line) and in haloes from the MARE NOSTRUM UNIVERSE 
(filled circles) for a sample of cluster-sized haloes in the mass 
range 6 - 7 X 1O 14 M /h. The shaded regions enclose 25% - 90% 
and 90% — 100% for a sample of M0KA haloes, while the three 
solid curves illustrate the same regions for the haloes in the MARE 
NOSTRUM UNIVERSE. 

3-4-2 Host halo concentration 

Recent studies combining strong and weak lensing of galaxy 
clusters have not only disagreed with predictions of the num- 
ber of arcs, but the observed clusters also seem to have Ein- 
stein radii larger than predicted in the ACDM cosmology 
(Zitrin et al. 2011a,c). This may emphasise that observed 
clusters are more concentrated than those found in numer- 
ical simulations (Oguri et al. 2005, 2009). However, over- 
concentrations in observed haloes may also result from an 
orientation bias (Meneghetti et al. 2011). 

We shall now quantify by how much the strong 
lensing signal changes if the normalization of the mass- 
concentration relation is increased. We perform different 
simulations of the same halo sample, assuming that the 
mass-concentration relation is given by Eq. (5), where we 
include a normalization factor Co. A reasonable choice for Co 
is a value between —3.2 and 6.4, consistently with simula- 
tion results for the mean concentration scatter at fixed halo 
mass. 

Figure 12 shows the median ratio of strong lensing cross 
sections with our fiducial sample as a function of c' vir — 
c V ir + Co for l/w — 5, 7.5 and 10. A higher normalization 
tends to increase strong lensing. Going from a mean value 
of the concentration for cluster size haloes from 4 to 8, the 
strong lensing cross section increases by about a factor of 4. 

3.5 Subhalo abundance 

Torri et al. (2004) have shown that merger events and sub- 
structures increase the strong lensing cross sections, mostly 
if the latter are located near the cluster centre. 

Our analysis has so far used cluster maps whose sub- 
halo mass function was obtained by sampling the distri- 
bution 7 down to 10 10 Mq/H. To statistically analyse how 
the strong lensing signal depends on the minimum subhalo 
mass in our models, we have generated a sample of triax- 
ial haloes with BCG, sampling the subhalo mass function 



Figure 11. Dependence of the median strong lensing cross sec- 
tion on the normalization of the mass-concentration relation. The 
solid line shows the median cross section for our fiducial sam- 
ple of cluster at 2; = 0.25, the dotted and the dashed lines 
show the relation for higher and lower value of the normaliza- 
tion in the mass-concentration relation relation: c' vir (M v i r , zi) = 
c v i r (M v i r , zi) + Co, for Zhao et al. (2009) cq = 0. 

down to a variable minimum mass m m i n . The total halo 
mass M V i r = 10 15 Mq//i is kept fixed. 

Figure 13 shows the strong lensing signal as a function of 
the minimum subhalo mass. The left panel shows the median 
ratio of runs with different resolutions compared to the fidu- 
cial case, whose subhalo mass is > 10 10 M©/ft. The strong 
lensing signal increases when the minimum subhalo mass 
changes from 1Q 10 M@/h to 10 14 Mq/H because the smooth 
halo mass tends to increase for larger value of m m i„, rais- 
ing the projected mass density distribution near the host 
halo centre. The right panel shows the ratio {a m m/o) at 
l/w = 7.5 as a function of the minimum subhalo mass ex- 
pressed in terms of the total mass. Again, the strong lens- 
ing signal increases because the smooth component contains 
more mass. This is shown in the same figure, where the data 
points show the median of the smooth mass component in 
the host halo as a function of the minimum subhalo mass, 
rescaled with respect to the smooth mass component, with 
a minimum subhalo mass of 10 10 M Q /h. 

Notice that the strong lensing cross section depends on 
the minimum subhalo mass. A decrease of the latter by a 
factor of 10 tends to increase (77.5 by 5% because the smooth 
mass increases. However, the impact of substructures can 
be different in different clusters depending of the particular 
configuration of the lens 

3.6 Generalized NFW density profile 

One of the main potential problems of the CDM model re- 
gards the inner slope of halo density profiles. Studying a 
sample of six strongly lensing clusters, Sand et al. (2004) 
concluded that, at 68 % confidence, the inner slope is con- 
sistent with p ~ 0.52 and that this is inconsistent with p = 1 
at the 99 % level. Recent analyses by Newman et al. (2009, 
2011) of Abell 611 and 383 have confirmed a flat central 
dark matter density profile. How much does the inner slope 
bias the strong lensing signal? 

To answer this question we have generated a sample of 



MOKA: Strong Lensing in Galaxy Clusters 11 




1 2 3 45678910 1 2 3 45678910 1 2 3 45678910 

C * + C C *ir + C C vir + C 



Figure 12. The median ratio of the strong lensing cross section with respect to the fiducial sample as a function of the normalization 
c v ir + co in the Zhao et al. (2009) mass concentration relation for three different value of l/w = 5, 7.5 and 10 in the left, central and 
right panels, respectively. The shaded region is enclosed by the lower and the upper quartiles. 




Figure 13. The dependence of the strong lensing cross section on the minimum subhalo mass. Left panel: median strong lensing cross 
section scaled by the fiducial triaxial simulation with subhalo mass resolution 10 10 Mq / h. Right panel: median rescaled cross section at 
l/w = 7.5 as a function of the minimum subhalo mass. The median rescaled strong lensing cross sections for l/w = 5 and 10 as a function 
of rrimin/Myir do not fall far away from the one estimated for l/w = 7.5. The data points connected with a dot-dashed line show the 
median of the host halo smooth mass component rescaled with respect to the smooth mass component when the minimum subhalo mass 
is 10 10 Mq /h, as a function of the minimum subhalo mass. 



haloes with different /3. We take adiabatic contraction in the 
center into account by means of Eq. (21), where mi(ri) is 
estimated by integrating equation (38). 
The generalized NFW density profile 



pgNFW = 



(r/r.)"(l+r/r.)3-^ ' 



(38) 



is taken to have an arbitrary inner slope (3. We define r_2 = 
(2 — j3)r s and the concentration c_2 = cw/(2 — /3) at the 
radius where the density profile is isothermal. 

Figure 14 shows the median strong lensing cross section 
for three samples of cluster sized haloes. NFW labels our 
fiducial sample with /3 = 1, while f3 = 1.5 and 0.5 illustrate 
a steeper and a shallower profile. 

The three panels in Fig. 15 show the relation between 
the median ratio of the strong lensing cross section at l/w = 
5, 7.5 and 10 compared to our fiducial sample as function of 
the inner slope /3. In each panel the shaded region encloses 
the lower and the upper quartiles, while the dashed line 
shows the least-squares fit to the data, 



log ((a(P)/a NFW )) = a ■ ft + b, 



(39) 
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Figure 14. The median value of the strong lensing cross section 
for three different values of the inner slope of the density profile 
of the host. Dotted, solid and dashed line refer to a sample of 
haloes with (3 = 1.5, 1 and 0.5, respectively. The thin solid curves 
show the quartiles of the NFW sample. 



for which we find 

a = 0.572 ± 0.005, 6 = 

a = 0.710 ± 0.005, 6 = 

a = 0.768 ± 0.009, 6 = 



-0.560 ± 0.005 a 5 
-0.701 ± 0.006 ar.s 
-0.750 ± 0.009 an, 
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A shallower (steeper) inner slope tends to decrease (increase) 
the strong lensing cross section with l/w > 5 by about 50% 
(100%). Consindering triaxial haloes with ft = 1 and /3 = 
1.5, analogous values for the strong lensing cross section for 
l/w — 7 and l/w = 10 and sources with z 3 < 1.25 have been 
found by Oguri et al. (2003). 



4 CONCLUSIONS 

We have presented a new algorithm to study the lensing sig- 
nal from cluster sized haloes. The ingredients that we use 
to build up our triaxial and substructured models take into 
account the most recent results from numerical simulations 
of structure formation. The possibility of turning these in- 
gredients on and off in our algorithm allows us to quantify 
the importance of all halo properties for their strong lensing 
efficiency. Starting from the halo deflection angle maps, we 
estimate strong lensing cross sections by ray-tracing. 
We can summarize our main results as follows: 

• Different structural halo properties affect strong lensing 
cross sections in different ways; 

• Averaging over a sample of cluster-sized haloes, we find 
that a central galaxy and triaxiality tend to enhance the 
strong lensing signal by 20% and 50 — 70%, respectively; 

• the strong lensing cross section monotonically increases 
with the host halo concentration; in a fixed mass bin, it is 
well characterized by log-normal distribution; 

• an increase (decrease) of the normalization of the mass 
concentration relation by a factor Co = 2 increases (reduces) 
the strong lensing cross section by a factor 5 (0.03); 

• increasing the minimum subhalo mass by a factor of 10 
weakly increases 07. 5 by about 3%; 

• the change in the inner slope of the density profile of 
the host halo is linearly related to y = log(o"i/u,); profiles 
shallower than NFW are weaker strong lenses; 

These results may help to understand future observational 
data and predictions for upcoming wide field surveys (Re- 
fregier et al. 2010). 
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APPENDIX A: SCATTER IN THE STRONG 
LENSING CROSS SECTION 

Haloes form hierarchically by gravitational instability of 
dark matter density fluctuations. Small systems form first 
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Figure Al. The median ratio of the strong lensing cross section 
between sample of haloes generated assuming a scatter ai n c and 
our fiducial one - where we have assumed o"i nc = 0.25. Different 
line types refer to different value of the variance in the log-normal 
scatter from which we draw the concentration to assign to each 
host halo in the respective sample. 



and merge to form larger objects. The assembly history de- 
pends on the environment where a halo grows. Virialized 
structures with the same final mass may have experienced 
different growth histories (Gao et al. 2005). Haloes with the 
same mass can thus have different concentrations and sub- 
halo mass functions. Numerical simulations have shown that 
the concentration distribution is well fit by a log-normal dis- 
tribution 



p(c\M v 



exp 



(lnc — In Cvir) 



2a. 2 

In c 



(Al) 



We assume <ri nc = 0.25 in our algorithm. In order to 
relate the concentration scatter with the scatter of strong 
lensing cross sections at fixed host halo mass, we have gen- 
erated different halo samples with different ai n c from which 
we randomly draw the host halo concentration. Figure Al 
shows the median ratio of the strong lensing cross sections 
compared to our fiducial sample for different standard devi- 
ations <7] nc . We find that smaller standard deviations reduce 
the median strong lensing signal. 

For each halo sample, we have estimated the stan- 
dard deviation <Ji og(T for length-to- width ratios l/w = 5, 
7.5 and 10. Figure A2 shows the relation between the stan- 
dard deviation in the strong lensing signal for haloes with 
M vir = 10 15 M Q //i, zi = 0.25 and z s = 2 and the scatter 
in concentration, rescaled to the scatter for the sample with 
fine = 0.25. For the chosen l/w ratios, the correlation is well 
fit by a straight line. Higher standard deviations in the con- 
centration scatter cause larger strong lensing signals. The 
dashed lines in all panels show the least-squares fit to the 
data points, 



(Tlogo- = a X (Tine + b, 

where we obtain the slopes 

a = 0.714 ± 0.049, for a 5 
a = 0.640 ±0.091, for a 7 . 5 
a = 0.506 ± 0.053, for a w . 



(A2) 
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Figure A2. The correlation between the scatter in the strong 
lensing cross section at three different value ofl/w and the scatter 
in concentration assumed to generate halo samples. 
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